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Abstract. At the New Horizons in Globular Cluster Astronomy conference (Padova, June 
2002), two members of the VLT globular cluster team presented different views on the im- 
portance of heavy-element sedimentation in Population II stars: "The lack of evidence for 
depletion of Fe and Li in the atmospheres of globular cluster subgiants led some people 
to susp ect that, for unkn own reasons, Population II stars are not affected by this mecha- 
nism." (Castellani 2003) and "There should be some mechanism that prevents sedimenta- 
tion." dGrattonll2003h . 

In this review, I will argue that the scepticism behind both these statements i s justified. We 
recent ly revisited the results on sedimentation in NGC 6 397 stars presented bv | Gratton et al] 
d200ll) using higher-quality VLT/FLAMES-UVES data dKorn et 311200651 2007). Element- 
specific abundance trends were identified which agree with atomic-diffusion predictions, 
turbulent mixing below the convective envelope is accounted for in a parametrized way. 
Have we thus detected signatures of Atomic Diffusion in Old Stars (ADiOS)? Or are these 
trends mere artefacts of Conspiring Inaccuracies in Abundance Observations (CIAO)? 
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1. On theory 

Since the pio neering work of 
lAller & Chapmanl {i960), microscopic 
diffusion has been studied in a variety of astro- 
physical applications, ranging from ordinary 
low-mass stars like the Sun to chemically- 
peculiar (CP) stars, horizontal-branch stars 
and white dwarfs. In the latter objects, effect 
of atomic diffusion (that is, sedimentation or 
levitation due to the net force atoms and ions 
will feel in stars) are large and well-studied 
(see Georges Michaud, these proceedings, 
and references therein). In solar-type stars, 
diffusive effects are considerably suppressed 
due to the presence of convection, but dif- 
fusive drainage of heavy elements from the 
convection zone can still occur. 



Early models (only treating gravitational 
settling, see Fig. 1) predicted sedimentation 
by a factor of 10 or more, in particular 
for me tal-poor stars with sh allow convection 
zones (Michau d et al.1 Il984l) . Apart from the 
potential impact on studies of Galactic chem- 
ical evolution, this prediction was relevant 
for the correct interpretation of the l ithium 
abund ances in warm halo stars dSpite & Spitd 
1982) which was identified with the amount of 
lithium produced in Big-Bang nucleosynthesis. 
Observational constraints on atomic diffusion 
in solar-typ e stars could only b e obtained in the 
mid-1990s dRyan et al.1 1 19961) and pointed to- 
wards near-negligible corrections for lithium. 

Over the years, the theoretical predictions 
were refined by including additional effects. 
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Fig. 1. Predictions of stellar-structure models including the effect of gravitational settling for a metal-poor 
([Fe/H]i n i, = -2) low-mass star that evolves to the turnoff in the course of 13.5 Gyr. A general depletion of 
the atmospheric abundances of heavy elements by factors of between 3 and 5 is predicted. 



Depending on the atomic structure, the down- 
ward diffusion of elements is more or less 
counteracted by the interaction of the atoms 
with the radiation field. This process, radia- 
tive acceleration, can levitate certain elements 
into the atmosphere, thereby raising their abun- 
dances above the original ones (see Fig. 2). 
Note that lithium is not levitated. This refine- 
ment was therefore incapable of addressing 
the observational fact that the Spite plateau 
of lithium is observed to be thin and flat, 
without a downturn towards hotter and more 
massive stars with shallower convection zones 
and larger predicted diffusion corrections. It 
was claimed that one simply had not ob- 
served a sufficient number of stars to de- 
tect t his predicted downturn ( Salaris & Weissl 
|2001|) . This idea did not really catch on 
with the observers, in spite of the accumu- 
lating evidence for a high baryon density £2b 
(and with it a high primordial lithium abun- 



dance) as inf erred from measureme nts of CMB 
anisotropics dNetterfield et al.ll2002l) . 

Around the same time, it became more 
and more evident that additional mixing pro- 
cesses are required to describe the properties 
of Am/Fm stars. In these, radiative accelera- 
tions accumulate iron at depths with temper- 
atures around 200 000 K leading to the cre- 
ation of an iron convection zone dRicher et al.l 
120001: iRichard et all |2001|) . To moderate the 
large abundance variations produced by unin- 
hibited atomic diffusion, one needs to intro- 
duce a mixing mechanism that keeps the lay- 
ers between the iron and the helium convection 
zone mixed. This process is referred to as tur- 
bulent mixing. It is parametrized by 

D T =400D He (:ro)(-^-r 3 
p(To) 

where Dj is the diffusion coefficient for tur- 
bulent mixing, DoeiTo) is the isotropic atomic 
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Fig. 2. Predictions of stellar-structure models including the effects of gravitational settling and radiative 
acceleration for the metal-poor star of Fig. 1. Depletions of up to a factor of 4 and enhancements of up to a 
factor of 2.5 are predicted by these models. Note the slightly shifted ordinate scale. 



diffusion coefficient of helium at the reference 
temperature Tq and p is the density. Thus, tur- 
bulent mixing is localized to a thin layer be- 
low the convection zone. This formulation was 
already applied to study the effects of tur- 
bulent jmxin^_pji_grayjtetional settling in the 
Sun dProffitt & Michaudl 1 1 99 l"h . For more in- 
sight into the physics of turbulent mixing, see 
Suzanne Talon (these proceedings, and refer- 
enc es therein). 

iRichard et ai](l2005l) showed that the inclu- 
sion of turbulent mixing allows to meet the ob- 
servational constraint of a flat and thin Spite 
plateau in the presence of significant surface 
depletions (see their Fig. 7): depending on the 
assumed strength (varying the reference tem- 
perature between log T = 6.0 and 6.28), deple- 
tions of 0.2 to 0.6 dex are predicted. 

In Fig. 3, such predictions are presented. 
As indicated above, turbulent mixing moder- 
ates both the large depletions and the enhance- 



ments. Some elements are now predicted to 
be depleted by * 0.2 dex (Mg, Al, Si), while 
others suffer a decrease of merely » 0.1 dex 
(Ca, Sc, Ti). Iron-group elements fall in be- 
tween these two (very moderate!) extremes. 
Identifying such abundance signatures is the 
challenge one is faced with. 

2. On observations 

It was already mentioned that the observations 
of the Spite plateau of lithium put some con- 
straints on uninhibited atomic diffusion. But 
atomic diffusion affects all elements to some 
degree, while the star evolves along the main 
sequence (MS). When the star finally becomes 
a red giant and the convection zone expands 
inward, the original composition is practically 
restored. Nature supplies laboratories where 
the effects of atomic diffusion can be stud- 
ied: globular clusters (GCs). One has to com- 
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Fig. 3. Predictions of stellar-structure models including the effects of gravitational settling, radiative accel- 
eration and T6.0 turbulent mixing for the metal-poor star of Figs. 1 and 2. Depletions by factors of 1.25 to 
2 are predicted by these models. Note the modified ordinate scale with respect to Fig. 1 and 2. 



pare atmospheric abundances of stars near the 
MS turnoff point (TOP) to those on the red- 
giant branch (RGB). Were systematic abun- 
dance differences in line with the theoretical 
predictions (like those of Fig. 3) identified? 

iKing et all (1 19981) used the Keck tele- 
scope to observe six subgiants in M 92 
([Fe/H] = — 2.3). While the resolving power 
was adequate (R = 45 000), the signal-to-noise 
ratio (S/N) was only 25-40. Because of this 
(and the lack of giant stars observed in the 
same program), they could not draw firm con- 
clusions about abundance differences between 
TOP and RGB stars. They state: "We note pos- 
sible evidence for [Fe/H] differences within 
M92T 

iGratton et ail d200ll) were the first to use 
UVES on the VLT to study unevolved stars 
in metal-poor globular clusters in greater de- 
tail. The main aim of their study was to ver- 
ify or refute the existence of elemental an- 



ticorrelations that had previously only been 
studied in GC giants. With S/N ratios around 
80 at fl = 40000, the nominal data quality for 
stars in NGC 6397 ([Fe/H] = -2) is good, but 
the spectra suffer from flat-fielding and order- 
merging artefacts which lead to systematically 
overestimated TOP-star e ffective temperatures 
(lKornetalJl2006bll2007h . As TOP and base- 
RGB stars were observed, the issue of atomic 
diffusion could also be addressed: "[Fe/H]'s 
obtained for TO-stars agree perfectly with that 
obtained for stars at the base of the RGB." 

ICohen & Melendezl d2005h studied 22 
RGB and three subgiant stars in M 13. They 
write in the abstract that "Most elements, 
including Fe, show no trend with r e ff [...] 
suggesting that [...] gravitationally induced 
heavy-element diffusion [is] not important 
[...]". However, a closer look at the pub- 
lished iron abundances reveals a systematic 
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Fig. 4. Chemical abundances (bullets and crosses) of FLAMES-UVES targets as a function of T e s- They 
are compared to predictions from models treating gravitational settling, radiative levitation and turbulent 
mixing, the latter with Tq varying between log T = 5.8 and 6.09. The intermediate T6.0 model described 
the element-specific trends remarkably well. 



difference between the two groups of stars: 
[Fe/H ]sGB = [Fe/H] R GB- 0.15. 

In iKorn et all (l2006al 120071) we presented 
the analysis of in total 18 stars covering the 
post-MS evolutionary sequence of NGC 6397 
in some detail: five TOP stars, two stars in the 
middle of the subgiant branch, five base-RGB 
stars and six RGB stars. They were observed 
with FLAMES-UVES on the VLT with nomi- 
nal data quali t y very similar to that achieved by 
lGrattonetal.ld20Qll) : R =48000, 70 < S/N < 
110. But since FLAMES-UVES is a fibre-fed 
spectrograph, the flat-fielding is more reliable 
and the intrinsic profiles of the strong TOP- 
star Ha lines can be recovered with confidence. 
In addition, we observed 130 stars along the 
subgiant branch with FLAMES-MEDUSA at 
somewhat lower resolution (R » 26 000). These 
spectra have a free spectral range of 200 A 



which facilitates Balmer-profile analyses sub- 
stantially. 

Note that it suffices to derive reliable 
stellar-parameter differences to study abun- 
dance differences between stars. For GC stars, 
log g differences can be accurately deter- 
mined from V-magnitude differences (with 
only weak dependencies on other parameters). 
Here, spectroscopic techniques are less pre- 
cise. For T e s differences, we rely on both pho- 
tometric and spectroscopic means. It is this 
combination of independent techniques that 
gives us confidence in our results. 

Based on both spectroscopic and photo- 
metric stellar parameters (see Table 1), we find 
trends of atmospheric abundances with evolu- 
tionary phase that are in good agreement with 
the predictions from models treating gravita- 
tional settling, radiative levitation and turbu- 
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Table 1. Stellar-parameter differences between TOP and RGB stars 



method A r eff (TOP - RGB) [K] A log g (TOP - RGB) A [Fe/H] 



Her 1124 ± 120 1.38 ±0.12 0.16 ±0.05 

V 1.38 ±0.05 

V - I 1070 ± 50 
B - V 906 ± 50 
v-y 1082 ±50 

b-y 1108 ±50 



lent mixing (see Fig. 4). The steepest trends are 
identified for Mg and Fe both of which are sig- 
nificant at the 3cr level. Similar trends are iden- 
tified among the 100+ FLAMES-MEDUSA 
stars using independent techniques (Lind et al., 
in prep.). 

3. On primordial lithium 

The observed lithium abundances of the TOP 
stars scatter around logs (Li) = 2.24 + 0.05 
(lcr). The two subgiant stars show a mean 
lithium abundance of logs (Li) = 2.36 + 0.06. 
Both the base-RGB and the RGB stars show 
lower lithium abundances indicative of dilution 
with lithium-free layers, as the convection zone 
expands inward (see Fig. 5). 

We can use the empirically calibrated 
atomic-diffusion model to predict the 
lithium abundances the NGC 6397 stars 
had at the time of formation. This value 
(log s (Li) p = 2.54 ± 0.1) is in good agree- 
ment with the predictions of Big-Bang 
nucleosynthesis at Q b as derived by WMAP 
dSpergel et all 120071) . In iKorn et all d2006al) . 
we called this "A probable stellar solution to 
the cosmological lithium discrepancy". 

4. On antagonism 

Bonifacio et al. I (120071) criticized our work in 
that it "relies heavily on the adopted tempera- 
ture scale, which is plausible, albeit inconsis- 
tent with the cluster photometry. An increase 
by only 100 K of the effective temperature as- 
signed to the TO stars would remove the abun- 
dance differences [...]". As Table 1 shows, 
our spectroscopic effective-temperature scale 



is not at all inconsistent with the cluster pho- 
tometry. Rather, it agrees within lcr error bars 
with three of the four colour indices studied 
and is conservative in the sense that it points to 
a larger T e ff difference yielding smaller abun- 
dance differences. Moreover, raising the TOP- 
star effective temperatures by 100K dimin- 
ishes the abundance trends by typically 0.05 
(Mg) to 0.07 dex (Li). While this may make the 
small abundance difference between lithium in 
TOP and subgiant stars (0.12 dex) insignifi- 
cant, it only removes one quarter of the iden- 
tified abundance trend for Mg. 

Criticism of more ge neral nature was raised 
bv lAsplund et all d2006l) . The detection of 6 Li 
in the spectra of several metal-poor stars puts 
constraints on atomic diffusion and mixing, as 
6 Li is even more fragile than 7 Li. Some of the 
atomic-diffusion models (those with the most 
efficient turbulent mixing compatible with the 
Spite plateau: T6.25 and T6.28) destroy 6 Li 
very efficiently lowering its abundance by up to 
factors of 100. Consequently, the original 6 Li 
abundance (observed to be « 5 % that of 7 Li) 
would be higher than the 7 Li abundance. Our 
observations point to a low-efficiency model 
(T6.0) that depletes roughly equal amounts of 
both isotopes (see Fig. 3). Together with pre- 
MS burning, this may aggravate the 6 Li prob- 
lem. Note, however, that the realit y of the 6 Li 
detec tions are now questioned dCayrel et alJ 
12007b . 

5. On diffusion in M 67 

Recently, IVandenBerg et all d2007l) studied 
the solar- metallicity, old (4Gyr) open cluster 
M 67 using up-to-date hydrostatic model at- 
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Fig. 5. Observed lithium abundances in the TOP and subgiant stars compared to trends predicted by the 
T6.0 model including the effects of gravitational settling, radiative acceleration and turbulent mixing. Note 
that there is a diffusion signature between the TOP and subgiant stars. The light-shaded area indicates the 
original lithium abundances of th ese stars, while the dark-shaded one represents the primordial lithium 
abundance as inferred by WMAP ( Sper gel et all2007h . 



mospheres as outer boundary conditions to the 
stellar-structure models. They found that the 
observed blueward hook at the TOP (indica- 
tive of a transition in mass (M tr ) from stars 
without a convective core to those with a con- 
vective core) is a sensitive function of the as- 
sumed composition and can only be repro- 
duced by models with a high solar metallic- 
ity sim ilar to that given by lGrevesse & Sauvall 
d 1 998b . Atomic diffusion may be able to lower 
Mtr sufficiently (by « 0.04 M ) to resolve this 
newly identified discrepancy. 

6. On the future 

We are only beginning to study the effects 
of atomic diffusion and mixing on atmo- 
spheric abundances of Pop II stars. More the- 
oretical and observational work is required 
to identify the physical process(es) respon- 
sible for moderating the effects of uninhib- 



ited atomic diffusion. The inclusion of turbu- 
lent mixing in a layer just below the convec- 
tive envelope is remarkably successful in de- 
scribing the observed abundance trends. The 
description applied is, however, in no way 
unique and physically satisfying. Mass loss 
may be ca pable of producing a similar m od- 
eration (see lVauclair & Charbonnej(ll995l) and 
Mathieu Vick, these proceedings). In the com- 
ing years, the ADiOSQ team hopes to put ad- 
ditional constraints on the dependence of tur- 
bulent mixing (and mass loss) on stellar pa- 



1 The ADiOS team currently consists of Martin 
Asplund, Paul Barklem, Lionel Bigot, Corinne 
Charbonnel, Remo Collet, Frank Grundahl, Bengt 
Gustafsson, Karin Lind, Lyudmila Mashonkina, 
Georges Michaud, Nikolai Piskunov, Olivier 
Richard, Suzanne Talon and Frederic Thevenin. 
The effort is led by the author. Seven of us attended 
this conference. 
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rameters, both observationally (by studying 
other globular clusters) and theoretically (by 
means of sophisticated (hydrodynamic) stellar- 
structure models). Encouraging results point 
to the importance of inter nal gravity waves 
(ICharbonnel & Talonll2005l) . 

Given our present ignorance about the ex- 
tent of atomic diffusion as a function of, 
e.g., metallicity, one should exercise caution 
when interpreting the abundance pattern of 
extreme Pop II s tars like HE 1327-2326 
(FrebeletalJ|2005). While among low-mass 
stars atomic diffusion produces largest effects 
at low metallicities, it is an omnipresent phe- 
nomenon which we should strive to include in 
XXI. Century modelling. 

Acknowledgements. I wish to thank all the scientific 
friends and colleagues who have contributed to the 
ongoing ADiOS effort. In particular, I would like 
to thanks Olivier Richard who supplied some of the 
figures for this article. 
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